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A non-equilibrium dynamic mechanism for the allosteric effect 
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Allosteric regulation is often viewed as thermodynamic in nature. However protein internal mo- 
tions during an enzymatic reaction cycle can be slow hopping processes over numerous potential 
barriers. We propose that regulating molecules may function by modifying the nonequilibrium 
protein dynamics. The theory predicts that an enzyme under the new mechanism has different 
temperature dependence, waiting time distribution of the turnover cycle, and dynamic fluctuation 
patterns with and without effector. Experimental tests of the theory are proposed. 

PACS numbers: Valid PACS appear here 
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A prominent property of enzymes (protein catalysts) 
is that their catalytic activities can be regulated. En- 
zymes that are allosteric have two or more binding sites. 
Effects of ligand (effector) binding or reaction on one site 
can propagate to another distant catalytic site and affect 
its activity. Understanding the allosteric mechanism(s) 
is an important topic in structural biology. Conventional 
models assume that effector binding modifies the equilib- 
rium confomational distribution of allosteric proteins 1] . 
Recent proposed "dynamic models" emphasize entropic 
(or the accessible configurational space) changes due to 
effector-binding induced modification of protein fiuctua- 
tion patterns [3, H, 0) S • Close examination reveals 
that these models actually share some basic ideas with 
the conventional models [7[ . In summary allosteric regu- 
lation is generally believed to be "fundamentally thermo- 
dynamic in nature" @. For later discussions, we charac- 
terize these existing models as being driven by "thermo- 
dynamic regulation". However, here we argue that the 
above thermodynamic description may be incomplete, 
and propose an alternative "nonequilibrium dynamic reg- 
ulation mechanism" . This idea is inspired by experi- 
mental and theoretical studies on dynamic disorder, the 
phenomena that the "rate constant" of a process is ac- 
tually a statistical function of time due to slow protein 
conformational motions [1, [13, [O, [H, [11] ■ 
In this work we specifically examine allosteric regulation 
of enzymatic reactions. Furthermore, we consider the 
case of positive regulation (i.e., effector binding results 
in higher activity) unless specified otherwise. The cat- 
alytic site being regulated can be described by a few slow 
global conformational modes (here we assume one for 
simplicity) and local conformational changes involving 
atomic rearrangement (here refering as the reaction coor- 
dinates). For barrier crossing processes, a system spends 
most of the time near the potential minima, and the ac- 
tual barrier-crossing time is transient. Therefore, one can 



reduce the potentials further to one-dimensional projec- 
tions along the conformational coordinate, and approx- 
imate transitions along the reaction coordinate by rate 
processes between the one-dimensional potential curves. 
Similar description has been used in other contexts (e.g., 
protein motor studies Protein dynamics is affected 

by substrate binding. A minimal model representing the 
states of a catalytic site is: E (empty), S (substrate 
bound), P (product bound). Fig. la illustrates an exam- 
ple used in this work. The protein states are described 
by potential curves along the conformational coordinate 
with localized transitions between them. For an enzy- 
matic cycle, a substrate molecule first binds onto the 
catalytic site (E^S), then forms a more compact com- 
plex from which a chemical reaction takes place (S— >P), 
and finally the product is released (P— >E). In the more 
familiar discrete kinetic form, the overall process can be 
represented &s E + S ^ ES ^ ES* ^ EP ^ E + P, 
with E, S and P refer to the enzyme, substrate, and 
product respectively. Notice that in general, the opti- 
mal conformational coordinates for reactant binding, the 
chemical reaction, and product release i nay not be the 
same (as also suggested experimentally [15'|), and some 
conformational motion is necessary during the cycle. Dy- 
namics of the reduced system can be described by a set 
of over-damped Langevin equations coupled to Markov 
transitions = + f 1(^)7 where x repre- 

sents the conformational coordinate, Ui{x) is the poten- 
tial of mean force at a given substrate binding state, Q is 
the drag coefficient, and fi is the random fluctuation force 
with the property < fi{t)fi{t') >= 2kBTC,i8{t - t'), with 
ks the Boltzmann's constant, T the temperature. Chem- 
ical transitions accompany motions along the conforma- 
tional coordinate with x-dependent transition rates. For 
simplicity we leave the more general description, the gen- 
eralized Langevin equation [ll,[il,[i3) for future studies. 
The dynamics can be equally described by a set of cou- 
pled Fokker-Planck equations (here we only consider the 
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FIG. 1: A minimal model for the catalytic site (a)Tlie free 
energy curves represent three distinct catalytic site binding 
states (E, S, P) projected onto the conformational coordinate 
X. The inlet illustrates that the smooth potentials are actu- 
ally coarse-grained over rugged potential surfaces. Effector 
binding at a distant site may modify the roughness of the 
potentials . (b) Two classes of reaction pathways along a 
multi-dimensional surface. Path 1 is along the long but low 
barrier minimum energy path. Path 2 is a short high barrier 
corner-cutting path. 
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Where Di — kBT/Ci the diffusion constant, kij the tran- 
sition matrix element, and pi the probability density to 
find the system at position x and state i. For simplic- 
ity we dropped the dependence of Ui and pi on x in the 
above equation and later discussions. This is a unified 
framework for describing allosteric regulation. The ex- 
isting models can be regarded as special cases with the 
conformational coordinate discretized The "thermo- 
dynamic regulation" models assume that effector bind- 
ing at a remote site can affect the dynamics at the 
catalytic site by modifying Ui. While differing in de- 
tails, these models assume a quasi-equilibrium distribu- 
tion along the conformational coordinate, thus a ther- 
modynamic treatment is appropriate Q. However, pro- 
tein conformational fluctuations can be very slow (e.g. 
from ms to minutes 0, [H, [iB|)i which is comparable 
or even slower than the enzyme turnover time. Conse- 
quently, protein fluctuations may not have fully accessed 
the conformational space, and thus not be in equilib- 
rium. Variation of the dynamic properties along the con- 
formational coordinate can have a dramatic effect on the 
apparent protein activity. Our recent theoretical anal- 
ysis showed that the observed slow protein conforma- 
tional dynamics can be explained by rugged protein po- 
tential surfaces^l?]. During relative motions between two 
protein parts, numerous noncovalent bonds(residue pairs 
with electrostatic, hydrophobic/hydrophilic, steric inter- 
actions, etc.) may form and break with associated local 
conformational changes. These processes are in general 
uncorrelated with each other, which result in rugged po- 
tentials uM. The relative motions are then character- 



ized by hopping over numerous potential barriers (re- 
fer to inlet of Fig. la). For a potential curve with 
random ruggedness, Zwanzig showed that the barrier- 
hopping process can be approximated by diffusion along 
a coarse-grained smooth potential with an effective diffu- 
sion constant D — Dp exp(— (e/fcsT)^), where Dq is the 
bare diffusion constant, and e is the potential roughness 
parameter ^T^. The reported value of e is 2-6 ksT [20j . 
With Dq ~ lO^^CTO^/s, D can be reduced to lA^/s with 
e ASksT. Thus internal diffusion can be a rate lim- 
iting step for enzymatic reactions and in principle can 
be regulated by allosteric effects (see Fig. 1). Further 
studies are necessary to clarify the atomic view of the 
proposed potential roughness regulation. The effective 
diffusion can be accelerated by inducing local conforma- 
tional changes and synchronizing the breaking and for- 
mation of the noncovalent bonds. It may be also related 
to the coupling mechanism between global and local vi- 
brational modes dicussed by Hawkins and McLeish [5]. 
In addition to being a theoretical possibility, the nonequi- 
librium regulation mechanism also has the following im- 
provements over the conventional regulation mechanisms. 
First, it is a more effective way to regulate the enzyme 
activity than the thermodynamic regulation (conforma- 
tional and entropic) mechanisms. To increase the activity 
by 10^° through an Arrhenius process, the activation "free 
energy" barrier needs to be lowered by 23 fesT. Similar 
amount of free energy change is needed for an equilibrium 
population shift mechanism over the reactant conforma- 
tional space. On the other hand, for an internal diffusion 
limited process, the reaction rate is linearly dependent 
on the effective diffusion constant. To increase the activ- 
ity by the same 10'^°, the lower bound of the roughness 
parameter only needs to be adjusted by ~ Sfc^r. 
Secondly, compared to the conformational change mech- 
anism, the nonequilibrium dynamic regulation mecha- 
nism has less requirements on the mechanical properties 
of the protein, similar to the proposed equilibrium dy- 
namic models. The distance between the two binding 
sites of an allosteric protein can be far (e.g., 15 nm for 
the bacterial chemotaxis receptor [2lj). Under the con- 
formational change mechanism, effective coupling of the 
two sites requires a faithful flnely tuned transmission of 
the mechanical strain due to ligand binding from one site 
to another one through a set of mechanical stress relay- 
ing network. These network residues must have mechani- 
cal properties distinctive from other residues to minimize 
energy dissipation to the surroundings. Otherwise, a sig- 
nificant portion of the effector binding energy would be 
wasted. In other words, coupling between these relaying 
residues and others should be minimized. By compar- 
ison, under the current nonequilibrium or the existing 
dynamics regulation mechanisms, the effect of effector 
binding can be highly nonlocal. Effector binding may 
affect the other site by finely regulating local structures 
far away from that site. By modifying the effective dif- 
fusion constant, these local modifications may affect the 
dynamics along the conformational coordinate x in our 
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TABLE I: 1 Model parameters. Here (E, S, P) refer to empty, 
substrate bound, and product bound, respectively. All en- 
ergy units in this table are in reduced units, ksT = 1. For 
simplicity,!) _B = Ds — Dp. The prefactors of the rate con- 
stants are chosen so that the maximum values of ksE, kps, 
and kss sxe approximately 10, 1.5, and 5, respectively. 

S^E P^E E^S P^S E^P S^P 
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formalism through couphng to global modes. The latters 
are composed of collective motions of residues vifithin the 
catalytic site and those far from it. The effect manifests 
itself through larger root-mean-square deviation as ob- 
served in NMR, x-ray crystallography, and in molecular 
dynamics simulations 0, M, HH • The dynamics of the col- 
lective motions should be examined as well. 

In our numerical calculations, the potentials are cho- 
sen to be harmonic potentials (see Fig. la), Ui — 
\n,(x - xoif + U° ,K, = (1, 0.5, 0.4), xo^ = (1, -0.5, 0), 
and — (0, —1, —3). To model transitions between dif- 
ferent states, we also model the transition state potentials 
by harmonic potentials, Ujj = ^{{x - xlj)/LijY + U^^. 
The transition rate from state j to i is given by fc^ {x) — 

/c°j exp[(J7j(a;) — J//^ (a;))/fcsT] with parameters given in 
Table 1. By solving Eq. [l] the enzyme turnover rate 
(which measures how many substrate molecules can be 
transformed into product by an enzyme molecule per unit 
time) is calculated by r = / (fc_Bppp — kpEPE)dx. Eq. [T] 
was solved with the algorithm developed by Wang etal. 
at given values of D and T [2^. The algorithm dis- 
cretizes the conformational coordinate, and transforms 
the partial differential equations into a jump process over 
many discrete states with their normalized populations 
p (defined as the probability density integrated over the 
discrete regions) described in the form Kp = 0. The 
composite K matrix contains transitions along both the 
conformational and reaction coordinates (see the original 
paper for details). Fig. 2a shows the calculated rela- 
tive enzyme turnover rate as a function of the internal 
diffusion constant. While the diffusion constant is not a 
rate-limiting parameter at high D values (as compared to 
the chemical transition rates), at smaller values of D the 
turnover rate depends on the diffusion constant linearly 
which is a signature of the existence of diffusion-limited 
steps. Fig. 2b shows the temperature dependence of the 
turnover rate. With high values of D, the exponential 
1 /T dependence mainly comes from the Arrhcnius depen- 
dence of the transition rates. However with small values 
of D, the turnover rate shows strong non-exponential de- 
pendence, since the effective diffusion constant D has a 
Gaussian dependence on 1/T. Experiments can test the 
predicted different temperature dependence of enzyme 
activity with and without the effector provided it is reg- 
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FIG. 2: Theoretical predictions, (a) Enzyme turnover rate 
as a function of the effective internal diffusion constant D. 
(b) Temperature dependence of the enzyme turnover rate. 
D = Do exp[— (e/fcsT)^], where e is the roughness parameter, 
and Do = 10^. Solid line:e = 0. dashed line: e = 4fcsr, the 
dotted line gives the exponential dependence if the Arrhenius 
rate equation is assumed. The temperature dependence of 
Do is neglected in this calculation, (c) Turnover waiting time 
distribution with D — 1 (solid line) and D — 0.001 (dashed 
line) . 

ulated by the current mechanism. Fig. 2c also shows 
the waiting time distribution between two consecutive 
turnover cycles calculated using the formula derived by 
Gopich and Szabo [11], 

P(t) = ltVexp( / (K - V)dt)Vp/(ltVp), (2) 
Jo 

where 1^ refers to vector contraction, K the above com- 
posite transition matrix, V the transition matrix with 
only product release transitions as the sole nonzero ele- 
ments, p the obtained steady-state population distribu- 
tion. A system with low D values shows non-exponential 
distribution due to dynamic disorder. At high D val- 
ues, effects of the dynamic disorder diminish and the dis- 
tribution is exponential. Therefore, we predict that an 
enzyme functioning under the new dynamic regulation 
mechanism shows larger dynamic disorder effects. This 
can be tested by measuring consecutive single enzyme 
turnover time distributions with and without the effec- 
tor, an extension of the work done by the Xic group p^ . 

For an enzymatic reaction under allosteric regulation, 
effector binding changes its reaction rate from ki to 
/c2- Let's define an effective free energy barrier change 
AAGt = fcsTl In(fc2/A:i)|, which must be due to the ef- 
fector binding energy. Is it necessary that the effector 
binding energy be no smaller than AAG^? In a related 
question, for a system described by a multi-dimensional 
potential surface shown in Fig. lb, is it possible that 
path 1 is dynamically comparable or even unfavorable 
than path 2? The answer to the latter is yes, provided 
that path 1 involves slow diffusion processes so that the 
average time for the transition along path 1 is even longer 
than path 2. A similar situation is discussed for tunneling 
pathways (e.g., proton-transfer reactions): the so-called 
corner-cutting lar ge c urvature tunneling and small cur- 
vature tunneling [24|. In this work we discussed that 
slow diffusion within a protein is physically possible due 
to the rugged potential surfaces. In his barrier cross- 
ing theory, Kramers derived dependence of the barrier 
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crossing rate on the barrier height, the drag coefficient, 
and other potential parameters of the system [25j . While 
most enzymatic reaction studies focus on barrier height 
changes, modification of the drag coefhcient can affect 
enzyme activity (see also ^26]). For proteins with rugged 
potentials, here we propose that the effective drag coeffi- 
cient and thus protein activity can be tuned over a broad 
range by modifying potential roughness. In this case, the 
effector binding energy need not to be less than AAG^. 
The dynamic mechanism discussed in this work is re- 
lated to other studies discussing protein dynamic prop- 
erties and allosteric regulation by considering the effect 
of rugged potential landscapes [4|). However, there is 
also a fundamental difference. The current model treats 
the enzymatic reaction as a nonequilibrium problem in 
general. In their NMR studies of dihydrofolate reduc- 
tase catalysis, Boehr et al. shows that the internal con- 
formational motion is the rate-limiting step [isj . This 
system may provide a nice test system for the proposed 
mechanism. In addition, slow conformational dynamics 
has been observed for allosteric proteins further sup- 
porting the validity of our formalism. Some experimen- 
tal observations supporting the existing dynamic models 
are also consistent with the current model. Kim et al. 
and Popovych et al. proposed that the allosteric signal 



is transmitted through dynamical rather than conforma- 
tional changes [111, [23]. We expect that the allosteric 
mechanism of a given protein has contribution from both 
thermodynamic regulation (the conventional conforma- 
tional change mechanism pj] and the newly proposed en- 
tropic effect [2,[llli,i), and nonequilibrium dynamic 
regulation proposed in this work. Different proteins may 
differ on which effect is dominant. 
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